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The specific heat, magnetic susceptibility and ESR signals of a Na-deficient vanadate Na^VgOs 
(x =1.00 - 0.90) were studied in the temperature range 0.07 - 10 K, well below the transition point 
to a spin-gap state. The contribution of defects provided by sodium vacancies to the specific heat 
was observed. It has a low temperature part which does not tend to zero till at least 0.3 K and a high 
temperature power-like tail appears above 2 K. Such dependence may correspond to the existence 
of local modes and correlations between defects in V-0 layers. The magnetic measurements and 
ESR data reveal S =1/2 degrees of freedom for the defects, with their effective number increasing 
in temperature and under magnetic field. The latter results in the nonsaturating magnetization at 
low temperature. No long-range magnetic ordering in the system of defects was found. A model for 
the defects based on electron jumps near vacancies is proposed to explain the observed effects. The 
concept of a frustrated two-dimensional correlated magnet induced by the defects is considered to 
be responsible for the absence of magnetic ordering. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.50.Ee, 75.45.+J, 76.30.-v 

I. INTRODUCTION 

The sodium vanadate NaV^Os is a unique compound 
with ID-magnetic structure formed within V-0 layers 
due to the special arrangement of vanadium- oxygen or- 
bitals. An intensive study of this compound began 
when a phase transition into a dimerized spin state was 
discoveredcl. According to the present model, vanadium- 
oxygen bonds give rise to a ladder structure with one 
electron (and one spin) periwig which should be shared 
between two vanadium ionsaia. Electrons are localized at 
rungs due to the Coulomb repulsion in accordance with 
the dielectric state of the crystal. In the high temperature 
phase the average charge of each vanadium ion is +4.5. 
The phase transition into a spin-gap state at T ==36 K 



was at first interpreted as a spin-Peierls transition . . Fur- 
ther theoretical and experimental investigationaJa point 
to a charge ordering phase transition that results in a 
zigzag charge distribution along the ladders, with a dou- 
bling of the lattice period and with alternating exchange 
interactions. The latter is responsible for the spin-gap 
opening according to Ref. ^|. The spin-gap causes an 
exponential freezing of the spin part of the, magnetic sus- 
ceptibility observed below the transitionE2l. 

There are two reasons why the problem of defects in the 
spin-gap crystals is of great interest. First, soliton-type 
spin clusters may be formed around imperfect spins or 



nonmagnetic impuritiestlhl3. These clusters are meso- 
scopic multispin objects with a microscopic total spin. 
The multispin nature of such clusters was confirmed by 
ESR experiments in the Ni-doped spin-Peierls compound 
CuGeOalij. Second, the defects in a spin-gap (and hence, 
nonmagnetic) matrix can induce the long range antifer- 
romagnetic ordering due to correlation of spins in neigh- 
boring clusters. Such an induced antiferromagnetic or- 
dering was observed in .the dimerized phase of the spin- 
Peierls crystal CuGeOata and in the Haldane compound 
PbNi 2 V 2 8 n. 

The substitution of magnetic ions by other magnetic 
or nonmagnetic ions in NaV^Os is not yet reported. The 
only reliable method to embed magnetic defects into this 
spin-gap crystal is the creation of sodium vacanciesEZI. 
Each vacancy causes a 3d-electron shared by two vana- 
dium ions to couple with an oxygen ion, thus making a 
next-to-defect rung of the vanadium ladder empty (i.e. 
nonmagnetic). This kind of defects should therefore be 
equivalent to the diamagnetic dilution or cutting the 
dimerized S =1/2 chains. Nevertheless, one should note 
that two ladder rungs situated next to the sodium va- 
cancy are electrically equivalent which allows the electron 
to occupy either site or to jump between the rungs. 

The present work is devoted to the study of low tem- 
perature thermal and magnetic properties of Na-deficient 
NaV2 05 in order to search for long-range magnetic or- 
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dering stimulated by defects. It is also interesting to 
compare the behavior of jumping magnetic defects in 
this kind of a dimerized spin matrix to the behavior 
of strongly localized defects in a spin-Peierls compound 
CuGeC>3. No transition to an ordered phase was found 
down to a temperature of 0.07 K, however unusual mag- 
netic and thermal properties due to the defects were ob- 
served and a model of jumping defects to describe these 
properties has been developed. The problem of delocal- 
ized magnetic defects is suggested for more detailed the- 
oretical investigation. Further experiments to clarify the 
question are proposed. 

II. EXPERIMENT 
A. Samples 

The sodium deficient single crystals were prepared by 
the techniques described in Ref. [l7|. Single crystals of 
the stoichiometric compound were embedded in a large 
quantity of the Na-deficient powder and were heated for 
one week. The resulting content of sodium in the crystals 
after heating was controlled with accuracy of about 1 % 
by comparing the results of x-ray and magnetic measure- 
ments to those obtained for reference powder samples. 
The content of sodium in the powder was determined 
from the molar ratios of the initial reagents. 

B. Specific heat 

The specific heat of Na x V205 samples with different 
sodium deficiencies was measured in the temperature 
range 0.3 - 8 K. The samples were put onto the holder 
mounted by thin thcrmoinsulating wires inside the mas- 
sive ring equipped by a thermometer. The temperature 
difference between the holder and the ring was controlled 
by Au-Fe thermocouple. The current induced by ther- 
moelectric power of the thermocouple was detected by 
SQUID. The signal from the SQUID was used in a feed 
back circuit to support the temperature difference equal 
to zero by heating the ring. Supplying the known power 
P to the holder with the sample we measured their tem- 
perature in real time, thus obtaining the total heat ca- 
pacity. More detailed description of the installation and 
the experimental technique including the correction for 
the parasitic power etc. is given in Ref. The specific 
heat of a stoichiometric sample was found to be cubic 
in temperature according to previous measurementslij. 
The molar heat capacities of the imperfect samples ap- 
peared to be much larger in the whole temperature in- 
terval. The contribution of the defects to the heat ca- 
pacity obtained as a difference between the total value of 
the molar specific heat of several Na^A^Os samples and 
that of x=1.00 sample is shown in Fig. 1. The specific 
heat of a perfect (x=1.00) sample fitted by T 3 -line is also 



given as a reference. The defects contribution may be di- 
vided into two main components: the low-temperature 
part remaining non-zero till at least 0.3 K (see Fig. 2) 
and the high-temperature power-like tail. The tempera- 
ture dependence above 4 K is approximately quadratic. 
The low temperature part increases monotonously in the 
concentration of defects (i.e. in 1 — x) while the value of 
the quadratic part reaches its maximum at 1 — x ~0.04. 
A test measurement of the specific heat under magnetic 
field was performed above 4 K. It was found to be field 
independent up to 7 T with an accuracy of 10%. Even 
at the largest vacancy concentration no sign of a phase 
transition into an ordered state was observed, in contrast 
to similar measurements on diluted CuGeC>3. 

The contribution of defects to the entropy was calcu- 

T 

lated numerically by taking the integral E = J A ^ F dT. 

0.3 

The temperature dependence of this integral divided by 
the molar entropy of an ensemble of two-level systems 
is shown in Fig. 3. One should note that the total en- 
tropies of the samples with x =0.96 and 0.92 are almost 
the same while the entropy of the x =0.98 sample is twice 
as small. This probably indicates the appearance of some 
collective state in the system of defects at high deficiency 
concentrations. 



C. Magnetic susceptibility 

In order to look at magnetic properties in the low tem- 
perature range, we performed a series of magnetization 
measurements using a SQUID magnetometer equipped 
with a miniature dilution refrigerator in the temperature 
range 0.07 - 1.2 K and in magnetic fields up to 8 T. 
The temperature dependencies of the magnetization at 
H ~0.05 T were obtained for samples with x=1.00, 0.98 
and 0.90. The corresponding inverse susceptibilities are 
shown in Fig. 4. The susceptibility of the stoichiomet- 
ric sample was found to be Curie-like, with the effec- 
tive concentration of free S =1/2 spins being equal to 
6T0 -4 . The small positive Weiss constant 9 ~0.03 K 
corresponds to weak antiferromagnetic interactions be- 
tween spins. The magnetic defects give rise to an addi- 
tional Curie-like contribution to the susceptibility which 
is nevertheless much smaller than that expected for the 
system of free 5 = 1/2 spins with the number equivalent 
to one half of the concentration of defects. The Weiss 
constant increases to 8 ~0.07 K. The effective concentra- 
tions of free S=l/2 spins and Curie- Weiss constants of 
several samples determined by fitting the magnetic sus- 
ceptibility below 1 K to a Curie law are given in Table 1. 

Table 1 



1 — x 9, mK Effective concentration 

"30±T0" 6-1Q- 4 

0.02 70±10 3.7-lQ- a 

0.1 70±10 6-10- 3 



2 



The inverse susceptibility for samples with small (1—x) 
values deviate from a linear temperature dependence for 
T >1 K. This seems to indicate that the effective concen- 
tration of spins increases with temperature. On the other 
hand, we cannot rule out that some of the temperature 
independent part may be due to non-uniformity along 
the sample holder (no corrections have been made to the 
present data). 

Test magnetization measurements were made for 
x =0.98 sample by standard Quantum Design SQUID 
magnetometer in the temperature range 4 - 77 K. The 
value of Curie constants obtained from data fits at 4 - 15 
K corresponds to (1 — x)/2 concentration of defects which 
is in agreement with previous resultsO. 

The field dependence of the magnetization at 
T =0.077 K was also studied for these samples. Although 
the temperature dependence of the susceptibility at low 
temperature is Curie like, the magnetization was found 
to have a component which does not saturate till H =8 T 
(see Fig. 5). It should be mentioned that at T =0.077 K 
the saturating field resulting in 90%-polarization of free 
5=1/2 spins is equal to 0.17 T. 

There was no indication of a phase transition down to 
0.077 K in accordance with the specific heat measure- 
ments described above. 



D. ESR 

The sodium deficiency gives rise to an increase in the 
ESR intensity at low temperature. Imperfect samples of 
NaV205 demonstrate the broadening of the ESR line in 
the temperature range below 3 K (see Fig. 6) in contrast 
to the constant value of the linewidth in stoichiomet- 
ric samplesEd. The ESR linewidth of magnetic defects 
caused by Na- vacancies appears to be narrower than that 
in a perfect sample which proves the different nature of 
the residual defects in a nominally stoichiometric sample 
and of the artificial defects in Na-deficient samples. 

III. DISCUSSION 

As mentioned in the introduction, the sodium defi- 
ciency in NaV2C>5 crystals appears to give rise to a new 
kind of magnetic defect in spin-gap systems. Unlike the 
usual nonmagnetic dilution, when impurity ions are fixed 
in the basic magnetic matrix and all the spins are local- 
ized, the absence of a sodium ion leads to the absence of 
an electron (but not an ion) at a rung of the vanadium 
ladder. Because the sodium atoms are positioned sym- 
metrically in between two rungs (see Fig. 7) one can see 
that the resulting unoccupied rung will have two almost 
equivalent positions. We suggest that this extra degree 
of freedom for the electron position is very important for 
the system with exchange interactions and we will try 
to explain the observed effects starting from this point. 



We also assume that the ground state of the electron is 
localized at one of the two next-to- vacancy rungs. In 
spite of the obvious loss in energy due to electron local- 
ization (of the order of the hopping amplitudeQ tn ~0.15 
eV), this assumption seems plausible. Mentioning, that 
the charge ordering phase transition observed at T=36 K 
is also followed by localization of electrons (at the ends 
of the rungs, with the corresponding hopping amplitude 
t± ^0.35 eV) one can suggest the existence of an addi- 
tional localization mechanism. For example, the shift of 
the electron to one of the rungs may break the symmetry 
of its potential hole in the vicinity of defect, thus strongly 
diminishing the hopping amplitude. 

The most obvious consequence of the sodium (elec- 
tron) deficiency may be understood neglecting for the 
exchange alternation. As mentioned above, the diamag- 
netic dilution of ID-magnets usually cuts magnetic chains 
into fixed segments with an equal probability for them to 
contain an odd or an even number of spins (we shall 
refer them to as "odd" and "even" segments). How- 
ever electron jumps can result in a disproportion between 
them. The ground state energy of a uniform even seg- 
ment is E even = —NEq — a/N and that of an odd one is 
E odd = -NE + b/N, where E ~ 0.886 J, b^ 2a ~ 3.5 J 
and N is the number of spins in a segmented. Thus, the 
even states correspond to a gain in exchange energy. A 
one electron jump between "odd-odd" and "even-even" 
states of two neighboring segments can change the en- 
ergy by approximately AE ~ (a + b)/N. Naturally, 
for long chains this effect should be strongly damped by 
dimerization which, nevertheless gives an additional gain 
in energy corresponding roughly to the spin gap energy, 
which comes from the recovering of " dimers" in even seg- 
ments. It is also reasonable to suggest the existence of 
a random electrostatic potential in a crystal with defects 
which makes the neighboring rungs slightly inequivalent. 
For simplicity one can assume it to be constant but of a 
random sign U — ±e. Taking the statistical distribution 
of segment length N, f(N) = (1 — x)x N (1 — x - con- 
centration of chains breaks) we obtain a set of two-level 
systems with the following molar heat capacity: 

C=\( 1 - x) n a E (°( A n) + C(A+)) , (1) 

where C(A^) is the heat capacity of a two-level system 
with the gap 

=| AE(N)±e |. This formula describes qualitatively 
the low temperature part of the experimental curves for 
various concentrations of defects (see Fig. 1,2). 

Consider also the behavior of such a system under mag- 
netic field which splits the energy levels corresponding to 
the odd states of the chains. One can easily obtain the 
following formula for the magnetization: 

tiff 

M(H,T) = -(l-x) f ,N A ^2f(N) ^ 6TT _~ 6 Z Em 

e~£T + e ~rt + e 

(2) 



3 



where the summation is performed over N for +e and 
— e. The result is in surprisingly good agreement with 
the corresponding experimental data (see Fig 4,5). This 
effect has a simple interpretation: due to a gain in ex- 
change energy, the number of even segments is initially 
larger than that of odd ones which is the reason for the 
system of defects to be almost nonmagnetic in the ground 
state. The effective number of magnetic objects (odd 
chains) may increase either by thermal activation or un- 
der magnetic field due to a gain in Zeeman energy. (Note, 
that the small steps at theoretical curves on Fig. 5 re- 
sult from the discrete distribution of two-level systems 
by energies AE(N)). One must emphasize that this very 
rough model does not cover the effects of dimerization 
and correlations between defects and ladders. In addi- 
tion it implies that the electrons are strongly localized 
at one of two next-to- vacancy rungs. The aim of this 
model is only to demonstrate the main features of the 
thermodynamic behavior of jumping defects. 

The observed low temperature widening of the ESR 
line in Na-deficient samples may result from the freezing 
out of electrons jumps in the vicinity of defects. When 
jumps are activated the ESR line should be narrowed by 
this kind of motion similar to the magnetic resonance line 
in a liquid. 

The unexpected result of the specific heat measure- 
ments is the observation of a power-like contribution at 
relatively high temperatures. While the low temperature 
part can be described in terms of the thermal activation 
of local modes (electron jumps) the power-like tail could 
indicate correlations of some kind between defects which 
may be of at least two origins. 

a) . The correlation between edge spins in segments due 
to a gain in exchange energy for "even" chains: 
The jump of one edge electron of a segment affects 
the electron at the other edge forcing it to shift in 
the same direction (for even number of rungs be- 
tween the vacancies) or in the opposite direction 
(for odd number) in order to recover the initial 
parity (see Fig. 7). Thus, the initial jump ap- 
pears to be a perturbation transmitting along the 
ladder. These shifts also disturb the zigzag pat- 
tern in neighboring ladders because of the correla- 
tions between them, which allows this perturbation 
to transmit also in the transverse direction. One 
should mention that Fig. 7 shows the perturba- 
tions damaging the dimerized state of the segments. 
It is also possible to imagine a situation where the 
dimerization is restored in which case it is obviously 
not the same gain in exchange energy. A detailed 
theoretical study of the correlations between such 
defects is required. 

b.) The spin clusters arising around the defect in a 
dimerized spin-gap matrix: Several spins should 
be correlated antiferromagnetically in the vicinity 
of each defect and the interchain interaction may 



cause the long-range AFM ordering)!!!! Such an 
ordering forces each pair of spins separated by an 
empty site to be parallel, otherwise the interchain 
interaction appears to be frustrated (see Fig 8). In 
case of NaV205 this frustration may be caused by 
AFM exchange J* between edge spins appearing 
due to electron jumps. It is probably the reason 
for the absence of the spin vacancy induced long 
range ordering confirmed in the present work down 
to T =0.077 K. Nevertheless, the defects should 
induce a frustrated 2D-magnet instead of an anti- 
ferromagnet. This is a strongly correlated system 
(see, for example, Ref. ^1]) whose specific heat is 
known to be power-like in temperature as observed 
in our experiment. 

It would be interesting to search for the phase tran- 
sition into an ordered state in crystals of NaV2 05 with 
a partial substitution of V-ions for other ions (e.g for 
Ti) fixed in the lattice. First, the susceptibility resulting 
from such a substitution would probably be much larger 
than those reported in this work. It may also result in the 
creation of either an ordered state or a frustrated mag- 
net, or a spin glass state. It will thus be very interesting 
to compare the properties of fixed and jumping defects 
and their influence on the spin-gap phase of NaV205. 
From the theoretical point of view it is very interesting 
to study the problem of localization of electrons in the 
vicinity of unoccupied rungs of the spin ladder (with and 
without dimerization) . The problem of the interaction of 
jumping electrons with the spin gap matrix has not been 
studied either. 



IV. CONCLUSIONS 

The contribution to the specific heat, magnetic suscep- 
tibility and ESR signal provided by defects in the spin 
gap matrix of NaV205 was observed. The data may be 
qualitatively described considering the defects as unoc- 
cupied rungs of a ladder structure. These defects pro- 
vide spin and jumping degrees of freedom. The temper- 
ature dependence of the specific heat reveals local modes 
and correlation between defects. We ascribe local modes 
to electron jumps in the vicinity of unoccupied rungs of 
the ladder structure. The correlation of defects is ex- 
plained qualitatively by a model of perturbations caused 
by electron jumps transferred along the ladder and to 
the neighboring ladders. Another possible mechanism 
of the correlation is the generation of magnetic clusters 
around defects which result in a frustrated 2D-magnet. 
The absence of an induced long range antiferromagnetic 
order (in contrast to other spin-gap systems) is inter- 
preted as the consequence of the geometric frustration 
of the exchange interaction between antiferromagnetic 
clusters. Magnetic properties of the system of defects 
are characterized by the temperature dependent effective 
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spin concentration and by the nonsaturating magnetiza- 
tion curves. These properties are also described in the 
frame of the jumping defect model. 
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Figure captions 

Fig. 1. The contribution to the specific heat from de- 
fects at various concentrations. Solid line is the quadratic 
high temperature approximation, dotted and dashed- 
dotted lines are calculations by formula (1) with the pa- 
rameters e=7 K, a + 6=450 K and 350 K for x=0.96 and 
0.98 respectively, dashed line is T 3 -fit to the phonon part 
of the specific heat (perfect sample) subtracted from all 
the data. 

Fig. 2. The residual low temperature part of the con- 
tribution to the specific heat from defects. Lines are the 
same as in Fig. 1. 

Fig. 3. The temperature dependence of the molar 
entropy of the defects. 

Fig. 4. The temperature dependence of the inverse 
susceptibility for samples with various Na-deficiency, H || 
6. Solid lines are theoretical fits by formula (2) with the 
parameters e=8 K, a + 6=400 K and 300 K for £=0.98 
and 0.90 respectively. Dashed line is the linear fit of low 
temperature points. 

Fig. 5. The field dependence of magnetization at 
T =0.077 K, H || 6. Solid lines are theoretical fits by 
formula (2) with the parameters as in Fig. 4. Dashed 
line is the magnetization of a paramagnet with the effec- 
tive spin concentration (1 — x)/2 (x=0.98). 

Fig. 6. The temperature dependence of the ESR 
linewidth at /=36 GHz, H || 6; solid lines are guide- 
to-eyes. 

Fig. 7. Schematic representation of the charge ordered 
ladder structure in the (a6)-plane of the NaV2C>5 crystal 
according to Ref. 0. Solid lines represent the paths of 
exchange interaction. Vanadium ions are positioned at 
the intersections of ladder rungs and legs; closed circles 
are V 4+ ions, open circles are V 5+ ions, grey circles are 
ions with jumping electrons; triangles are the projections 
of sodium vacancies onto the (a6)-plane. The pairs within 
dashed ovals are coupled by larger exchange interaction 
(" dimers" ) . The arrows show the direction of correlated 
electron jumps excited at finite temperature. 

Fig. 8. An antiferromagnetic interchain frustration 
due to electron jumps (frustrated region is outlined by 
dashed rectangle). For simplicity the zig-zag patterns 
are shown as chains. 
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